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ABSTRACT 

Aims. We attempt to determine the molecular composition of disks around young low-mass stars in the p Oph region and to compare 
our results with a similar study performed in the Taurus-Auriga region. 

Methods. We used the IRAM 30 m telescope to perform a sensitive search for CN N=2-l in 29 T Tauri stars located in the p 
Oph and upper Scorpius regions. *^CO J=2-l is observed simultaneously to provide an indication of the level of confusion with the 
surrounding molecular cloud. The bandpass also contains two transitions of ortho-Fl 2 CO, one of SO, and the C'^O J=2-l line, which 
provides complementary information on the nature of the emission. 

Results. Contamination by molecular cloud in '^CO and even C'^O is ubiquitous. The CN detection rate appears to be lower than for 
the Taurus region, with only four sources being detected (three are attributable to disks). Ft 2 CO emission is found more frequently, but 
appears in general to be due to the surrounding cloud. The weaker emission than in Taurus may suggest that the average disk size in 
the p Oph region is smaller than in the Taurus cloud. Chemical modeling shows that the somewhat higher expected disk temperatures 
in p Oph play a direct role in decreasing the CN abundance. Warmer dust temperatures contribute to convert CN into less volatile 
forms. 

Conclusions. In such a young region, CN is no longer a simple, sensitive tracer of disks, and observations with other tracers and at 
high enough resolution with ALMA are required to probe the gas disk population. 
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1. Introduction 


Star formation occurs in clusters that are embedded in collaps¬ 
ing molecular clouds (Lada & Lada 2003| l. The different stel¬ 
lar populations and properties, such as age, mass or density, can 
have effects on the circumstellar disk evolution, which strongly 
depends on the environment. In regions of high stellar density, 
disks can be more easily photoevaporated by the intense UV ra¬ 
diation Pohnstone et al.|1998| l. This means that the evolution of 
disks in the vicinity of massive stars will differ from the evo¬ 
lution of isolated disks (see, e.g., Mann et ^|2014 1. The disk 
structure and chemistry can therefore vary from region to region. 

Circumstellar disks are the sites where planets form and 
their gas and dust provide the raw materials for planet build¬ 
ing. Constraining the chemical evolution of protoplanetary disks 
therefore is a major challenge in understanding the planet for¬ 
mation process. Molecular-line emission is an important tool for 
deriving the disk characteristics ( |Dartois et al. 2003 [ Qi et al. 
|2006| l. Because H 2 , the most abundant molecule, cannot be ob¬ 
served, we rely on lower abundance tracers, principally simple 
molecules, to derive those characteristics. In a previous study 
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of ~ 45 stars in the Taurus region, we showed that CN N=2- 
1 is a good tracer of disks because it is readily detectable in at 
least 50 % of the observed sources ( |Guilloteau et al.|2013| here¬ 
after Paper I). Contamination by molecular clouds was minimal 
in this line; it was found in fewer than 10 % of the cases, but 
some sources exhibited strong lines that most likely originate 
from outflows rather than circumstellar disks. 

However, the Taurus-Auriga area is a region of isolated star 
formation, and the previous result may not be applicable for 
other star-forming regions, either because of stellar density or 
more simply, because of age. We extend here our study to the p 
Oph star-forming region, which is both younger and denser than 
the Taurus region. We report a search for CN, ortho-H 2 CO, SO, 
'^CO, and C'^O emissions in 22 young stars located in the p 
Oph region and 7 stars in the upper Scorpius area. We attempt to 
compare these results with previous observations performed in 
the Taurus region (Paper I), and more particularly, the CN con¬ 
tent in these star-forming regions. 

The paper is organized as follows. The observations and data 
analysis are described in Sect. 2. We present the results and the 
source analysis in Sect. 3. In Sect. 4 we discuss the difference in 
molecular composition observed between the Taurus region and 
that of p Oph. We conclude in Sect. 5. An Appendix displays all 
the observed spectra. 
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Object 

a 

(J2000) 

6 

(J2000) 

5y(3.4) 

(mJy) 

Sy(1.3) 

(mJy) 

Other 

Name 

Sp Type 

L. 

(Le) 

Mass 

(Me) 

Age 

(Myr) 

rms 

(mK) 

V1146SCO 

15:57:20.0 

-23:38:50.0 




MO 




28 

J1603-1751 

16:03:23.7 

-17:51:42.3 




M2 




36 

AS 205A 

16:11:31.4 

-18:38:26.0 

27.2 

450 

HBC 254 

K5 

4.0 



67 

SR 4 

16:25:56.1 

-24:20:48.3 

4.4 

31 

AS 206 

K5 

2.17 

1.14 

1.1 

41 

GSS 26 (*) 

16:26:10.3 

-24:20:54.9 

24.2 

125 


MO 

1.39 

0.56 

0.5 

34 

EL 20 

16:26:18.9 

-24:28:20.2 

7.3 

50 

VSSG 1 

MO 

0.93 

0.62 

1.1 

54 

LFAM 1 

16:26:21.7 

-24:22:50.8 

17.5 

250 






37 

DoAr 24E 

16:26:23.4 

-24:21:00.7 

8.3 

70 

GSS 31 





35 

DoAr 25 

16:26:23.7 

-24:43:14.1 

25.0 

280 

WSB 29 

K5 

1.43 

1.12 

2.1 

40 

EL 24 

16:26:24.1 

-24:16:14.0 

48.8 

390 

WSB 31 

K6 

2.58 

0.96 

0.6 

38 

EL 27 

16:26:45.0 

-24:23:08.2 

38.7 

300 

GSS 39 

MO 

0.78 

0.58 

1.2 

38 

WL 18 

16:26:49.0 

-24:38:25.7 

3.1 

85 

GY 129 

K7 

0.3 



70 

SR24S 

16:26:58.5 

-24:45:37.1 

26.6 

530 

HBC 262 

K2 

4.4 



78 

SR 21 

16:27:10.2 

-24:19:12.9 

4.2 

95 

EL 30 

G3 

11.38 

1.97 

2.2 

33 

J1627-2451B 

16:27:15.1 

-24:51:38.8 




M2 




28 

IRS 41 

16:27:19.3 

-24:28:44.4 

6.2 

< 60 

WL3 

K7 

1.61 

0.80 

0.8 

50 

CRBR 85 

16:27:24.7 

-24:41:03.2 

1.5 

150 



1.4 



22 

YLW 16c 

16:27:26.5 

-24:39:23.4 

6.5 

60 

GY 262 

Ml 

1.11 

0.48 

0.5 

23 

IRS 49 

16:27:38.3 

-24:36:58.8 

4.4 

25 

GY 308 

MO 

1.02 

0.64 

1.0 

47 

DoAr 33 

16:27:39.0 

-23:58:19.1 

3.7 

40 

WSB 53 

K4 

1.81 

1.44 

2.9 

38 

WSB 52 

16:27:39.5 

-24:39:15.9 

10.2 

51 

GY 314 

K5 

0.95 

1.04 

4.2 

41 

IRS 51 

16:27:39.8 

-24:43:15.0 

12.7 

110 

GY 315 





83 

Elying Saucer)*) 

16:28:13.7 

-24:31:39.1 

3 



Ml 

0.14 



40 

WSB 60 

16:28:16.5 

-24:36:58.0 

15.3 

89 

YLW 58 

M4 

0.23 

0.20 

0.9 

31 

SR 13 

16:28:45.3 

-24:28:19.2 

10.0 

60 

HBC 266 





37 

DoAr 44 

16:31:33.5 

-24:27:37.7 

10.4 

105 

HBC 268 

K3 

1.55 

1.29 

5.1 

38 

RNO 90 

16:34:09.2 

-15:48:16.9 

7.6 

25 

HBC 649 

G5 

10.24 

1.87 

2.3 

31 

Wa Oph 6 

16:48:45.6 

-14:16:36.0 

10.3 

130 

HBC 653 

K6 

2.32 

0.98 

0.7 

26 

AS 209 (*) 

16:49:15.3 

-14:22:08.6 

17.5 

300 

HBC 270 

K5 

2.11 

1.18 

1.2 

24 

HD 163296 (*) 

17:56:21.3 

-21:57:22.0 




Al 


2.25 

5.0 

50 


Table 1. Source sample. Sources in upper Scorpius are listed in italics. Sources marked with (*) are the sources for which CN 
emission from disk has been detected. The epoch 2000 coordinates are adopted from SIMBAD. Sy is the continuum flux at 3.4 mm 
and 1.3 mm. Spectral types and luminosities (given for ~ 130 pc) are taken from Ricci et al. p010|l ex cept for AS 205A (Bast et al. 
201 l|l, WL 18 ([Andrews et al.|2010| , SR 24S ( Andrews et al.|2010| l, CRBR 8 5 ([Pontoppidan et al.||2005| l, and the Flying Saucer 


( [Grosso et al.||2003| ). Masses and ages are derived from the tracks reported in jPalla & Stahler| ( |1999| l using the luminosities and 
effective temperatures from Ricci et al. (2010 1 . 


Table 2. Detected CN emission 


Source 

Line Flux 

Velocity 

Width 

name 

(Jy kms^*) 

(kms^*) 

(kms“‘) 

AS 209 

2.1 ±0.3 

4.5 ± 0.4 

4.2 ± 0.7 

GSS 26 

2.2 ± 0.3 

4.0 ± 0.2 

2.2 ± 0.4 

IRS 41 

1.5 ±0.3 

5.5 ±0.1 

0.9 ± 0.2 

CRBR 85 

^ 2.5 

see 

FigH 

Elying Saucer 

3.4 ±0.5 

4.1 ±0.3 

4T ± 0.6 


Notes. CN line parameters for the detected sources. 


made when we selected the sources. Sources with outflows were 
deliberately included. The distribution of sources is displayed in 
Fig{^ which shows the location of the studied sources overlaid 
on an extinction map. Of these 30 stars, 22 are located in the 




Oph star-forming region, for which a reliable distance, 120 pc 
based on VLBI astrometry ( Loinard et al.|20()8a] l or optical ex 
tinction (Knude & Hog|1998|l, is available. The distances to the 


other 8 stars might lie in the range 120 to 170 pc (Loinard et al. 
'2008b de Geus et al.|198^: we assume 160 pc hereafter. We re 


port our observations of these stars, but cannot reliably analyze 
them in terms of stellar luminosity and mass. 


2. Observations and data analysis 

2.1. Source sample 

Table [T] lists the 30 stars involved in this study. HD 163296, an 
isolated HAe star, was used for calibration purposes. Our main 
objective was to search for potential targets for the stellar mass 
determination, to test models of early stellar evolution. Hence, 
our sample is biased toward late spectral types. 

The overall sample covers a wide range of spectral types 
(from M4 to Al), luminosities (from 0.2 to 12 Lq), and 1.3 mm 
continuum flux densities (from < 4 to ~ 800 mJy), but is not 
complete with respect to these quantities. No specific account 
of extinction and location with respect to molecular clouds was 


2.2. Observations 


Observations were carried out with the IRAM 30 m radiotele¬ 
scope on October 12-13, 2012, and from August 8-14, 2013. The 
weather conditions were quite varied, with water vapor content 
changing from 3 mm on good nights to 6-7 mm on others. Given 
the low elevation of the sources (mostly below 25°), this resulted 
in single-sideband system temperatures ranging from 350 K to 
about 1000 K, with averaged values between 400 K and 600 K 
for most sources. However, pointing conditions were good. 


However. 


Loinard et al. 


1 2008b I reported for the two stars WL 5 and 


VSSG 14 (located near IRS 41) preliminary measurements indicating a 
distance of ~ 165 pc. 
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50m00s 40m00s 30m00s 20m00s lOmOOs 16h00m00s 



30s 28m00s 30s 16h27m00s 30s 26m00s 

Right Ascension (J2000) 


Fig. 1. Extinction maps of the p -Oph star-forming region (based 
on the 2MASS PSC |Dobashi et al.|2013| l, contours from 1 to 18 
mag in steps of 1 mag. The bottom map is a zoom on the area 
in which most sources are located (black square in the top map). 
Asterisks are sources with CN detection. 


A few strong lines of SO, SO 2 , and SO^ found in one of 
the sources allowed us to check the sideband rejection, which 
was found to be around 13 dB. All observations were made at 
relatively low but constant elevation. We thus did not apply any 
gain-elevation curve correction and used a simple uniform con¬ 
version factor from antenna temperature (T^ to flux density of 
Jk - 9 Jy/K. A comparison with previous observations on HD 
163296 (Paper I) suggests that our calibration is accurate within 
10 %. 

The observations were performed in symmetric wobbler¬ 
switching mode, with the two references +60" away at the same 
elevation. The weather conditions were stable enough to provide 
flat baselines. 

As the system temperature varied much from day to day, we 
observed the sample of 30 sources (see Table [T]) for integration 
times varying between 1 and 8 hours. The resulting sensitivity at 


the nominal resolution of 6V ^ 0.26 km.s * is given in Table 
(rms column). For a typical line width of AV = 3 km.s *, the cor¬ 
responding Icr error on the integrated line flux can be obtained 
by multiplying this number by Jk - VAVdy ^ 8 Jy.km/s/K. 
Our survey depth is thus ranging between 0.18 to 0.65 Jy km.s"\ 
with a typical value of 0.3 Jykm.s^'. 

2.3. Data analysis 

We applied the same analysis procedure as in Paper I, remov¬ 
ing linear baselines in a window + 60 km.s ' wide around each 
spectral line, and fitting Gaussian profiles to detect lines, with 
proper account for the hyperfine structure of CN and C'^O. 

The two transitions of ortho-H 2 CO differ only by the small K 
splitting, with the 211.2 GHz line on average 20% stronger than 
the 225.2 GHz line. We computed an average H 2 CO spectrum 
and line intensity using this ratio. 


3. Results 


As described in Paper I, line emission from these objects can 
come a priori from four distinct regions: molecular cloud(s) 
along the line of sight, the circumstellar disk, a molecular out¬ 
flow, and (especially for the youngest objects) a remnant en¬ 
velope. A molecular cloud will exhibit narrow lines, and be¬ 
cause of our observing technique, which differenciates the emis¬ 
sion from that +1' away, may appear in emission or absorp¬ 
tion. Disks should lead to symmetric, in general double-peaked, 
line profiles, and line widths oc sini (with typical values 2-4 
km.s *, face-on objects being rare). Contribution from outflows, 
envelopes, or molecular clouds may vary considerably and may 
be difficult to distinguish from disk emission, so that our results 
can be dominated by contamination. 


Figure|^shows the spectra toward the two illustrative sources 
CRBR 85 and the Flying Saucer, an edge-on disk at the periphery 
of p Oph ( [Grosso et al.|2003 1. Spectra for all other sources are 
presented in the Appendix. 

The '^CO results are completely dominated by contamina¬ 
tion. Only four sources show no obvious sign of contamination: 
AS 205, J1603-1751, DoAr 44, and VI146 Sco. No disk com¬ 
ponent can be fit to the observed profiles, except for the already 
known case of the isolated HAe star HD 163296. 


Even in C'^O J=2-l, contamination exists, with negative sig¬ 
nals in IRS 41, IRS 49, El 20, SR 21, El 24, El 27, and WSB 60, 
and strong emission toward CRBR 85, LFAM 1, YW 16c, and 
GSS 26. We have a weak detection in HD 1632926 (see below). 

On the other hand, CN emission is extremely rare. The only 
emissions that can be attributed to disks are from HD 163296, 
AS 209, the Flying Saucer (see Fig|^, and perhaps GSS 26. Our 
line flux for AS 209 is slightly lower tha n that reported from 
SMA observations by Oberg et al.| ( |2011 1 , probably as a result 
of contamination, and we also detect the second group of hyper- 
fine components that was below the sensitivity threshold for the 
SMA. We have clear evidence of cloud contamination in DoAR 
24E, LFAM 1, because of negative features, and in CRBR 85. 
In the latter case (Fig|^, the spectrum is the superposition of a 
very narrow component, with highly anomalous hyperfine ratios, 
and a broader component that may be due to an underlying disk. 
The case of IRS 41 remains unclear: there is some CN emission 
(1.1 + 0.3 Jy km.s^'), but relatively narrow, ~ 0.9km.s“*, and at 
a velocity where C'^O shows residual contamination. Line pa¬ 
rameters for the detected sources are indicated in Table |2] 
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'^C0(2-1) 

C'^0(2-1) 

CN(2-1) 

HyCO 

50(54-43) 


Fig. 2. Lines toward CRBR 85 and the Flying Saucer. Top panel; '^CO J=2-l spectrum. Bottom panel, from top to bottom: '^CO 
J=2-l, C'^O J=2-l, CN N=2-l, the average of both 0 -H 2 CO transitions, and SO 54 - 4^, shifted in intensity by multiples of 0.15 K 
to avoid overlap. The intensity scale is the antenna temperature (T^ in K): conversion to flux density can be obtained using a factor 
of 9 Jy/K. The spectral resolution is 0.26 km.s''. For the Flying Saucer, the bottom panel spectra has been smoothed by a factor 2 
for better clarity. 


For H 2 CO, contamination also seems to be rather 
widespread: negative signals are observed in El 20, SR 21, SR 
4, El 27, DoAr 24E, and TEAM 1. Positive signals are detected 
in GSS 26, WSB 60, El 24, YEW 16c, and CRBR 85. This is a 
ratio of 11 to 29 sources. 

SO is rare; there is strong emission in TEAM 1, and a strong 
negative signal in IRS 41, but nothing elsewhere. 


3.1. Specific objects 

HD 163296. Although only included here as a calibration 
source, we note that we have a 5cr detection of C'^O, with a total 
line flux of 4.3 +0.9 Jy km.s'*. The CN and '^CO line fluxes are 
consistent with previous measurements, indicating good calibra¬ 
tion. 


GSS 26. This is one of the possible new disk detections in CN. 
The line width, however, is very narrow, 0.75 ± 0.12km.s'*, 
and is in fact consistent with that derived from H 2 CO lines, 
0.63+0.06 km.s'*, although the systemic velocity differs slightly 
(3.81 + 0.06 vs 3.43 + 0.03km.s''). Either this is a disk seen 
nearly face-on, or the detected lines are residual contamination 
by the surroundings, as may be suggested by the rather strong 
C'^O emission. 


CRBR 85. The spectra of CRBR 85 (also known as 
CRBR2422.8-3423, see Comeron et al.|1993t Pontoppidan et al. 


2(105 1), clearly show two components: a very narrow compo- 
0.5km.s'*) superimposed on a broader emission 
• 2km.s'*). The narrow component exhibits very 


nent (AV ~ 

(AV ~ 1.5 

unusual hyperflne ratios for CN N=2-l, which, if interpreted as 
saturation, would indicate very high opacities. However, they are 
most likely the result of our differential observation technique 
in a cloud extending over 2 arcmin that is nearly opaque in the 
strongest hyperflne components, as for CW Tau in Paper I. 

The broad component may be attributable to disk or enve¬ 
lope emission. As CRBR 85 is a highly inclined, embedded ob¬ 
ject, the limited linewidth indicates a very low mass star if this 
emission comes from a disk. The other possibility is emission 
from a larger distance, that is, from the envelope. The low line 
flux would then indicate low excitation, since the detection of 
weak hyperflne components indicates substantial opacity. Note, 
however, that the main group of hyperflne components may be 
seriously affected by absorption or contamination due to the nar¬ 
row component, which would result in higher apparent opacities 
than in reality. 


LFAM1. Also known as GSS 30 IRS 3, this is a strong radio 
source with emission at 6 cm ( |Leous et al.||1991|). This source 
is classified as a Class I source (Bontemps et al. 2001) l. The 
spectra in this direction are very unusual. We have strong emis¬ 
sion from SO and SO 2 , but also the ion SO^, as well as sev- 
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eral isotopologues of SO and SO 2 . H 2 CO is quite bright, and 
CN shows a negative feature that indicates contamination by a 
molecular cloud. S-bearing molecules are in general enhanced 
in shocks. In fact, SO^ is a clear indicator of shock chemistry 
( Turner||1992|l, predict ed to be abundant in dissociative shocks 
( Neufeld & Dalgarno||1989) l, and has recently been detected in 
the jet-driven outflow of LI 157 ( |Podio et al.|2014[ ). The unusual 
spectrum of LFAM 1 therefore probably indicates an outflow 
from a very young object. Another possibility is some contri¬ 
bution from an outflow driven by GSS 30 IRS 1 because LFAM 
1 is located 15" to the northeast of this Class I source, on the 


edge of its extended V-shaped nebula (see Fig. 8 of Allen et al. 
2002|l, typical of a large outflow cavity. 


YLW 16c. This source is the only one clearly showing the 
double-peaked spectrum expected from Keplerian rotation in 
one molecule, here in H 2 CO. However, the H 2 CO/CN ratio is 
high. In Paper I, we only had strong H 2 CO emission and/or high 
H 2 CO/CN ratio from sources with outflows, or relatively mas¬ 
sive stars, while YLW 16c is of spectral type Ml. 


4. Discussion 

4.1. p Ophiuchi and Taurus: two different regions 

The properties of sources in p Oph look very different from those 
in Taurus-Auriga: only two potential disk detections in CN (GSS 
26, Flying Saucer), plus perhaps CRBR85 and one source in 
SO (LFAM 1) out of 22 sources. In Taurus, there are 21 disks 
detected in CN (this number excludes the sources that also have 
SO emission of unknown origin) out of 46 sources, and 8 sources 
emitting in SO out of 41. 

Our survey depth ranges between 0.18 to 0.65 Jykm.s"*, 
with a median value of 0.32 Jy km.s *, a mean of 0.36, and a dis¬ 
persion of 0.17. Rescaled to the same distance (i.e., corrected by 
a factor (120/140)^), the Taurus observations of Paper I reached 
a median sensitivity of 0.20 Jykm.s^* (mean 0.25) with a dis¬ 
persion of 0.17. There is thus a typical factor of 1.5 in sensitivity. 

However, this does not appear to be the cause of the large 
apparent difference in detection statistics between the two re¬ 
gions. Figurej^shows the CN line intensity vs. the 1.3 mm con¬ 
tinuum flux (rescaled to a common distance of 120 pc) for both 
regions. Red indicates that CN has been detected. The Taurus re¬ 
gion exhibited a high disk detection rate (> 50 %), whereas only 
three or four sources have been detected in the p Oph region. 
Furthermore, the population of bright disks found in Taurus has 
not been detected in our p Oph study. 

We also performed a more detailed study of the CN flux 
density distribution in both regions. We applied a statistical 
treatment to our data that takes upper limits into account. We 
used the Kaplan-Meier estimator, which provides a maximum 
likelihood distribution for a randomly censored dataset. A part 
of our data being limited (upper limits), this estimator allows us 
to estimate the fraction of sources for which the flux is below (or 
above) a given value. The details about this method can be found 
in [Feigelson & Nelson] ( |1985| l. This estimator only requires 
that the censoring he random, which is the case in our samples 
since we have considered sources lying at different distances. 
We computed the Kaplan-Meier estimator for detections and 
nondetections (upper limits) for both samples. The results are 
displayed in Fig|^ The distributions show the probability that 
a given source in our sample has a CN flux higher than the 



Fig. 3. CN line intensities vs 1.3 mm continuum flux. Red dots 
highlight detected sources in p Oph and Taurus, while blue is 
for sources in upper Scorpius. All fluxes have been rescaled to a 
distance of 120 pc. 


abscissa value. The difference between the two samples is large: 
for the Taurus region, the median value for the flux is ~ 0.8 
Jykm.s^*, whereas this value is about twice lower for the p 
Oph sample. Only 20 % of the sources have a flux higher than 
1 Jykm.s^* in p Oph, whereas the sources in Taurus are much 
brighter since 40 % of the sources have a flux higher than 1 
Jykm.s ’. 


4.2. Possibie cause of these differences 

The CN disk detection rate clearly is much lower in the p Oph 
region than in the Taurus region. Furthermore, we observed a 
lack of bright sources in p Oph. We discuss in this section some 
hypotheses that might explain the differences between the two 
regions. 

Different stellar properties (sample bias)? As discussed in 
Sect 2.1 the two samples were not built using the same criteria. 
Unlike the Taurus sample, which included a variety of stars from 
spectral type M4 to AO, our p Oph sample mostly includes stars 
in the M2-K3 range and has no warm star (only one G5 and one 
G3). However, CN is most easily detected precisely in this M2- 
K3 range in Taurus sources (Paper I). If we restrict the Taurus 
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Fig. 4. Cumulative flux density distribution for CN emission in 
Taurus (blue line) and p Oph (red line), constructed using the 
Kaplan-Meier estimator to include upper limits and scaled to a 
common distance (d=120 pc). The sample size are 46 stars for 
Taurus (Paper I) and 22 stars for p Oph. 


sample to these spectral types, the difference between p Oph and 
Taurus is even more striking. 


poorly known. Given all uncertainties, it is possible that viscous 
spreading in disks located in the Taurus region produces disk 
sizes larger by 1.2 to 2 times than those in p Oph. Even a limited 
size difference might be sufficient to explain the observed dis¬ 
tribution; according to Fig|^ 80 % of the sample sources have a 
flux below 1 Jy km.s * in p Oph, while for Taurus this is obtained 
for 2-3 Jykm.s^*, which means that radii differing by 1.4 - 1.7 
could bring the two distributions to agree reasonably well, p Oph 
would still lack the brighter sources found in Taurus, but these 
are presumably outflows, not disks. 

The smaller size might also be a result of initial conditions 
linked to the more crowded environment in p Oph, an effect that 
known to affect the size of the protostellar envelopes (|Motte 


et al.|1998 1 . The missing protostellar material at large distances 


is the one with the highest angular momentum, which means that 
disks may be born smaller in p Oph. 


Different CN chemistry? The intrinsic differences in disk 
properties between p Oph and Taurus can affect the CN chem¬ 
istry and invalidate the assumption of similar CN surface densi¬ 
ties, however. Because p Oph is younger, the stars are brighter; 
from evolutionary tracks, a factor of about 2-3 in luminosities for 
stars between 1 and 3 Myr is expected. This in turn should lead 
to a factor (2 - 3)'^^^ = 1.2 - 1.3 in disk temperature at the same 
radius. Moreover, disks are not at chemical equilibrium, there¬ 
fore age matters. In the next section, we present the results we 
obtained with our chemical model that evaluates the importance 
of these effects. 


Contamination? Might contamination be the only cause of 
the difference? While it is obviously stronger in this region, this 
looks very unlikely. Contamination by clouds is best traced by 
H 2 CO and C'^O, which always show very narrow lines. Such 
narrow lines would not have prevented a detection of the wider 
lines coming from Keplerian disks in CN. 


Smaiier disks? In Paper I, we suggested that CN surface den¬ 
sities were essentially constant (except in the warmer sources). 
Under this assumption, which was substantiated by chemical 
modelling ( Chapillon et al.|[2012| l, the CN line flux provides a 
nearly direct measure of the disk size. The apparent radii derived 
from the line flux under this simple hypothesis agreed well with 
interferometric determination in Taurus. With the same assump¬ 
tions, the measured flux densities indicate here that we have de¬ 
tected three disks with outer radii 500 - 600 au (GSS 26, AS 209, 
and the Flying Saucer), have 3cr upper limits around 500 au for 
three other sources (WF 18, SR 24S and IRS 51), and that most 
sources are smaller than 320 au, with only one putative disk at 
450 au, IRS 41. 

The smaller size observed here may be due to the source 
ages, disks in p Oph being younger and therefore having less 
time to spread by viscous effects. Making this statement quanti¬ 
tative is difficult, however, because the age difference between p 
Oph (age < 1 Myr |Fuhman & Rieke||1999|l and Taurus (~ 1 
Myr, but with sources up to a few Myr Fuhman et al. 2003| l 
is not well known. In the framework of self-similar viscously 
spreading disks ( |Fynden-Bell & Pringle||1974| l, the disk radius 
evolves with time as where y is the radial exponent of 

the viscosity and T - \ + t^lts, t, being the disk age and 
the viscous timescale. The classical a disk prescription yields 
y = 312 - q , where q is the radial exponent of the temperature 
law, thus y ~ 1 as ^ ~ 0.5 for typical disks (see, e.g., |Guilloteau| 
|et al.|20lT] for derivations). Viscous timescales are unfortunately 


4.3. Impact of the chemistry 


We used the Nautilus gas-grain model described in Semenov 
|et al.| ( [20T0| l and |Reboussin et JI] ( |2014[ ). Nautilus computes the 
abundance of species as a function of time in the gas phase and 
at the surface of the grains. The chemical network contains 8624 
reactions; 6844 are pure gas-phase reactions 1780 are grain- 
surface and gas-grain interactions. The gas-phase network used 
for this work is the same as described in |Reboussin et al.| ( |2014| l. 
The full network is available on the KIDA (Kinetic Database 
for Astrochemistry) website The dust grains are represented 
by spherical particles with a radius of 0.1 pm and are made of 
amorphous olivine. The gas and dust temperatures were assumed 
to be the same, and we used a dust-to-gas mass ratio of 0.01. To 
obtain the initial chemical composition of the disk, we first com¬ 
puted the chemical composition of the parent molecular cloud. 
We ran Nautilus during 10^ yr for typical dense cloud conditions; 
a gas density of 2x10"^ cm“^, a temperature of 10 K, a visual 
extinction of 10, a cosmic ray ionization rate of 1.3x10“'^ s"', 
and a C/O ratio of 0.5. The elemental abundances used for this 
work are listed in [Reboussin et al^ ( |2014] l (see Table 1 of their 
paper), except for the oxygen elemental abundance, for which 
we considered the value 3.3x10“^ (a low-depletion case). For 
the disk physical parameters, we used the parametric disk model 
described in jHersant et al. (20091 except for the vertical temper¬ 
ature profile, which is based on Williams & Best p014| l. The 


model corresponds to a disk mass of 0.03 Mq with a l/r‘-^ sur¬ 
face density profile out to r = 700 au. The disk mid-plane tem¬ 
perature was used as a free parameter. The abundances of chem¬ 
ical species in the disk were calculated up to 5x10® yr. 

Figure |5 shows the CN column densities as a function of 
the disk midplane temperature at a radius of 300 au, at three 


http://kida.obs.u-bordeaux 1 .fr/models 
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Disk Midplane Temperature (K) 


Fig. 5. Panel (a) represents the two extreme temperature profiles at 300 au for a weak temperature gradient, and the corresponding 
CN column density as a function of the disk midplane temperature for three different ages; 1, 2, and 5 Myr are shown in panel (b). 
Panel (c,d): as (a,b), but for a strong temperature gradient. 


different disk ages and for two different temperature profiles 
(small or strong vertical gradient). Overall, the CN column 
density decreases as the temperature increases and even reaches 
values well below any detectable levels (down to a few lO" 
cm“^) at temperatures above 40 K in the disk midplane. With 
the 30 m radiotelescope, our detection limit is around 3x10'^ 
cm“^ for a disk of 400 au radius for a detection level of 1 Jy 
km.s“'. This value does not vary much with the temperature. 
The temperature dependence is complex and non-monotonic 
because the dominant reaction routes for CN change with 
temperature. The gas-phase CN abundance is mostly affected by 
grain-surface chemistry: at higher temperatures, the diffusion 
of the species at the surface of the grain is much more efficient, 
which contributes to convert molecules on grains into less 
volatile forms. As this conversion process is relatively slow, the 
temperature dependence is more pronounced at later ages. The 
difference between the two temperature profiles is quite small at 
high temperatures, but for very low temperatures (below 15 K), 
the CN column density is lower when we use a strong tempera¬ 
ture gradient. These results agree with the observations obtained 
in Paper I, in which CN was undetectable in most warm sources. 


As mentioned before, all these simulations were performed con¬ 
sidering the same temperature for the gas and dust. Models in 
which these two quantities differ indicate that the decisive factor 
for CN chemistry is the dust temperature. Figure also shows 
that the 20-30 % higher temperature of disks in p Oph is suffi¬ 
cient to affect the CN surface density by a factor of a few, except 
for very low temperatures (below 10 K). Similar differences 
exist as a function of disk ages for disk temperatures above 15 K. 


5. Conclusions 

We have performed a sensitive survey of 29 young stars in the p 
Oph and upper Scorpius regions in CN, ortho-H 2 CO, SO, ^^CO, 
and C^^O rotational lines near 206 - 228 GHz with the IRAM 
30 m telescope. Compared to a similar study performed for the 
isolated star formation region Taurus-Auriga, the detection rate 
of CN is much lower in the p Oph region. 

This result may indicate that disks in the p Oph region are on 
average smaller than those in the Taurus-Auriga complex, per¬ 
haps because they have not spread out sufficiently by viscous 
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diffusion, or as a result of initial truncation of their parental pro- 
tostellar condensation. 

However, the CN chemistry is shown to be sensitive to disk 
temperature, so that a direct comparison of the disk proper¬ 
ties would require resolving the disks out instead of relying on 
the good correlation between CN flux and disk radii previously 
found in the Taurus region. Interferometric observations are re¬ 
quired for this purpose. 

Finally, these results show that even large single-dish tele¬ 
scopes are severely limited in identifying disks around embed¬ 
ded young stars because of contamination, but also because of 
sensitivity. Disks smaller than about 300 au, which represent 50 
% of the disks in the sample studied by Guilloteau et al. ( 2013| l 
in the Taurus region, are beyond the sensitivity limit of even the 
largest telescope operating at 1.3 mm, the IRAM 30 m. Only 
ALMA is sensitive enough in the southern hemisphere to detect 
a substantial fraction of the gas disk population. 
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Appendix A: Spectra for individual sources 

This appendix displays the spectra toward the various sources 
(in the same order as listed in Table [T]i. For each source, the 
top panel shows the '^CO J=2-l spectrum. The bottom panel 
displays on a common scale from top to bottom the spectra of 
'^CO J=2-l (with fit as in top panel), C'^O J=2-l, CN N=2-l, 
the average of both 0 -H 2 CO transitions, and SO 54 - 43 , shifted in 
intensity by multiples of 0.15 K to avoid overlap. The intensity 
scale is the antenna temperature (T^ in K): a conversion to flux 
density can be obtained using a factor of 9 Jy/K. 

The spectral resolution is 0.26 km.s^^ in the upper panel, 
while the other spectra have been smoothed by a factor 1, 2, 
or 4 for better clarity. 
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Fig.A.l. Lines toward V1146 Sco. 



Fig.A.l. Lines toward J1603-175L 
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Fig. A.3. Lines toward AS 205A. 



Fig.A.4. Lines toward SR 4. 
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Fig.A.6. Lines toward EL 20. 
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Fig.A.7. Lines toward LFAM 1 . 



Fig.A.8. Lines toward DoAr 24E. 
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Fig.A.9. Lines toward DoAr 25. 
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Fig.A.lO. Lines toward EL 24. 
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Fig.A.ll. Lines toward EL 27. 



Fig. A.12. Lines toward WL 18. 
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Fig. A.13. Lines toward SR 24S. 
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Fig. A.15. Lines toward J1627-245IB. 



Fig.A.14. Lines toward SR 21. 



Fig. A.16. Lines toward IRS 41. 
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Fig. A.17. Lines toward CRBR 85. 



Fig. A.18. Lines toward YLW 16c. 




Fig. A.19. Lines toward IRS 49. 



Fig.A.20. Lines toward DoAr 33. 
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Fig. A.21. Lines toward WSB 52. 



Fig. A.22. Lines toward IRS 51. 



Fig. A.23. Lines toward Flying Saucer. 



Fig. A.24. Lines toward WSB 60. 
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Fig. A.25. Lines toward SR 13. 



Fig.A.26. Lines toward DoAr 44. 
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Fig. A.27. Lines toward RNO 90. 



Fig. A.28. Lines toward Wa Oph 6. 
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Fig. A.29. Lines toward AS 209. 



Fig. A.30. Lines toward HD 163296. 
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